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ABSTRACT 

Classical Be stars are known to occasionally transition from having a gaseous circumstellar disk ( "Be 
phase" ) to a state in which all observational evidence for the presence of these disks disappears ( "normal 
B-star phase" ) . We present one of the most comprehensive spectropolarimetric views to date of such a 
transition for two Be stars, tt Aquarii and 60 Cygni. 60 Cyg's disk loss episode was characterized by a 
generally monotonia decrease in emission strength over a time-scale of ^^1000 days from the maximum 
V-band polarization to the minimum Ha equivalent width, consistent with the viscous time-scale of 
the disk, assuming a ^0.14. tt Aqr's disk loss was episodic in nature and occurred over a time-scale 
of ~2440 days. An observed time lag between the behavior of the polarization and Ha in both stars 
indicates the disk clearing proceeded in an "inside-out" manner. We determine the position angle of the 
intrinsic polarization to be 166.7±0.1° for tt Aqr and 107.7±0.4° for 60 Cyg, and model the wavelength 
dependence of the observed polarization during the quiescent diskless phase of each star to determine the 
interstellar polarization along the line of sight. Minor outbursts observed during the quiescent phase of 
each star shared similar lifetimes as those previously reported for fi Cen, suggesting that the outbursts 
represent the injection and subsequent viscous dissipation of individual blobs of material into the inner 
circumstellar environments of these stars. We also observe deviations from the mean intrinsic polarization 
position angle during polarization outbursts in each star, indicating deviations from axisymmetry. We 
propose that these deviations might be indicative of the injection (and subsequent circularization) of 
new blobs into the inner disk, either in the plane of the bulk of the disk material or in a slightly inclined 
(non-coplanar) orbit. 

Subject headings: circumstellar matter — stars; individual (pi Aquarii) — stars: individual (60 Cygni) 



1. INTRODUCTION 



Since the pioneering work of lStruvel (|l93lh . a sub-set of 
B-type main sequence stars rapidly rotati ng at velocities 



ranging from >40% (for early-type B stars: ICranmei 



2005 ) 



20041) 



to possibly 100% of the critical rate (jTownsend et al 
have been identified. Photometric, spectroscopic, polari 
metric, and interferometric observations of these so called 
classical Be stars indicate that they are characterized by 
the presence of a geometrically flattene d, circumstellar de- 
cretio n disk (for a recent review see iPorter fc Riviniui 
120031) . While for single stars the origin of the gas in 
these disks is agreed to be material ejected from the stel- 
lar photosphere, the precise mechanism for forming these 
disks remains elusive. Suggested disk formation mecha- 
nisms include the compression of s tellar winds ("WCD" 
model; iBiorkman fc Cassinellilll993l). potentially aided by 
the channeling of magnetic fields ( Cassinelli et al.l 120021 : 
iBrown et al.l 20081). as well as the viscous decretion disk 
model ( Lee et aTl 1991), which cou ld be fed via non-radial 
pulsational events~ ([Cranmeij|2009t ) similar to those which 
have been observed in several Be stars (jRivinius et al.l 
[19981 iNeiner et al.l[200l) . 

Classical Be stars are known to exhibit variability 
a multitude of time-scales and amplitudes, includ- 



ing short-term (hours to weeks) polarimetric variabil- 
ity (jCarciofi et al.l [20071 ) . multi-periodic short-term spec- 
troscopic variability arising from non-radial pulsations 
(jRivinius et al.l [19981 ) . and long-term (several to tens of 
years) cyclical spectroscopic variability o riginating frorn 
one-armed dens ity waves in the disk (jOkazakil 119911 : 
IStefl et al] l2009f) . Arguably the most dramatic type of 
variability observed is the aperiodic "Be to normal B to 
Be" transition whereby classical Be stars dissipate (i.e. 
lose) their disks, and will often begin to regenerate a 
new disk after a period of relative quiescence (see e.g. 
lUnderhill fc Doazanlll983 and references therein). While 
several individual stars have been observed to experience 
this transition one or more times (Underbill fc Doazan 



Il982t iDoazan et a"Llll983l : IClark et al.„200a : .Vinicius et al 



it is not known whether there is a characteristic 
time-scale for disk-loss and regeneration that would sug- 
gest a singular disk formation mechanism or whether each 
event is unique. Although their sam ple size is small, recent 
monitoring of eight open clusters bv lMcSwain et al.l (j2009f ) 
has revealed a surprising 12 new transient Be stars out 
of a sample of 296 stars monitored, suggesting that disk- 
loss/renewal episodes may be a prevalent phenomenon. Be 
stars experiencing a "Be to normal B to Be" transition are 
clearly of significant astrophysical importance as they pro- 
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vide observational evidence of how Be disks form, thereby 
constraining theoretical efforts to identify the disk forma- 
tion mechanism. 

Polarimetry has an established history of providing a 
unique diagnostic in identifying and studying the de- 
tailed circums tellar environments of both individual clas- 
sical Be stars (iQuirrenbach et allll997l : IWood et al.lll997t 
IClarke fc Biorkman' '1998*1 and larger statistical ensem- 
bles (ICovne fc Kruszewski 1969; McLean fc Brown 19?! 
iPoeckert et al.lll979t I Wisniewski et all l2007bl ) . Electron 
scattering modified by pre- or post-scattering absorption 
by hydrogen atoms in the disk produces the characteristic 
wavelength depe ndent polarization signature of Be stars 
(jWood et al.|[l996a ih). Modehng such data yields infor- 
mation about the density in the inner r egion of these disks 
(|Wood et al.lll997tlCarciofi et al.ll2009D and has been used, 
along with interferometric obse rvations, to confirm tha t 
the disks are geometrically thin (|Quirrenbach et anil997f ). 

The two stars we examine in this paper, 60 Cygni and 
TT Aquarii, are well known classical Be stars. 60 Cyg (HD 
200310) is a BlVe classical Be star (I SlettebakllT98ai . Pe- 
riodic radial velocity variations of Ha and He I 6678 A 
suggest that the star might be a spectroscopic binary hav- 
ing a period of 146.6 ±0.6 days, while 1.064 day line pro- 
file variations and 0.2997 day photometric variations may 
be caused by n on-radial pulsations ( Koubskv et al.ll2000f) . 
iKoubskv et al.l (|2000, ) present evidence that the star has 
experienced disk-loss and subsequent regeneration in the 
past (see e.g. their Figure 1), with photometric variability 
positively correlated with Ha emission strength; however, 
they neither discuss the detailed time-scale for these events 
nor publish their data in sufficient format to derive this 
informat ion, tt Aqr ( HD 212571) is a Bl Ill-IVe classical 
Be sta r (jSlettebaklll982ll . As discussed in iBiorkman et al.l 
(|2002[ ) and references therein, tt Aqr's most recent disk 
likely developed in the early 1950s and persisted until 1996. 
Analysis of spectroscopic data obtained after 1996, when 
the star was in a diskless "normal B star" phase, revealed 
the star has a clos e binary co mpanion with an orbital pe- 
riod of 84.1 days (jBiorkman e t al. 2002), which might be 
engaged in variable mass transfer with the primary. 

In this paper, we present and analyze an extensive spec- 
tropolarimetric dataset which covers the disk-loss phase 
transition in two Be stars. In Section [21 we outline our 
new spectroscopic and spectropolarimetric observations of 
60 Cyg and tt Aqr. We analyze the time evolution of our 
spectroscopic data in Section [3. II and the time evolution of 
our polarimetric data in Section 13.21 A thorough discus- 
sion of the interstellar polarization along the line of sight 
to 60 Cyg and n Aqr is presented in Section r3.3l Finally we 
offer a detailed discussion and interpretation of our results 
in Section 0] and suggest future work in Section [51 

2. OBSERVATIONS AND DATA REDUCTION 

The spectropolarimetric data analyzed in this study 
were obtained by the University of Wisconsin's (UW) 
HPOL spectropolarimeter, mounted on UWs 0.9m Pine 
Bluff Observatory (PBO) telescope, tt Aquarii was ob- 
served on 127 nights between 1989 August 8 and 2004 Oc- 
tober 10 while 60 Cygni was observed on 35 nights between 

^ IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the National Science Foundation. 



1992 August 3 and 2004 September 26 (see Tables[T]and[2l). 
Data obtained before 1995 were recorded using a dual Reti- 
con array detector spanning the wavelength range of 3200 - 
7600A with a spectral resolution of 25A (Wolff ct al. 1991). 
Since 1995, HPOL has been equipped with a 400 x 1200 
pixel CCD and two gratings which provide coverage from 
3200 -6020A and 5980 - 10,500A at a resolution of 7A be- 
low 6OOOA and lOA above this point (jNordsieck fc Harris! 
[ ^96) . Further deta i ls abo ut H POL can be found in Noolg 
(I1990D . imiff et all (|T99l . and lHarries eTall ((2000l) . and 
all HPOL data obtained between 1989-1994 are publically 
available from the Multimission Archive at STScI (MAST) 
website : http : / /archive . st sci . edu /hpol 

HPOL data were reduced, processed, and analyzed 
using REDUCE, a spectropolarimetric software pack- 
age devBlo2edbj;_t he University of Wisconsin-Madison 
(se dWolff et al.|[T996i ). Routine monitoring of unpolarized 
standard stars at PBO has enabled the instrumental po- 
larization to be carefully calibrated; the residual instru- 
mental systematic errors depend mildy on the date of the 
observations, but range from 0.027-0.095% in the U-band, 
0.005-0.020% in the V-band, and 0.007-0.022% in the I- 
band. HPOL spectroscopic data are not calibrated to an 
absolute flux level due to the nonph otometric skies typ i- 
cally present during the observations ([Harries et al.ll2000l ). 

We supplement our HPOL spectroscopic observations 
of 60 Cyg with 65 nights of spectra obtained from 1998 
August 26 to 2006 December 9 using the fiber-fed echelle 
spectrograph of the 1.0m Ritter Observatory at the Uni- 
versity of Toledo (Table [1]). These spectra span nine non- 
adjacent 70 A wide orders in the range 5285-6595A with 
R ~26,000. The data were reduced in IRAF using stan- 
dard spectroscopic techniques, including bias correction, 
flat-fielding, and wavelength calibr ation. A dditional de- 
tails about Ritter data are given in Morrison e t al. ( 1997t). 

Our HPOL polarimetric observations of tt Aqr are also 
supplemented b y literature V-band polarization m easure- 
ments (Table m iMcDa^ [19861 119901 Il99l . [T999I) . More- 
over, both our 60 Cyg and tt Aqr data are supplemented by 
unpublished observations made with the Lyot Polarime- 
ter at PBO betwen 1979-1987 (Table [3]), which spanned 
an effective wavelength range of 4800-7600 A in 1979, 
4800-8200 A in 1980, and 4600-7200 A in 1987. Fur- 
th er details about the PBO Lyot Pol arimeter can be foun d 
in lLupie fc Nordsieckl (|T987ll : Whitnev fc ClavtonI (|1989[ ): 



and [Nook et al] ([1990[ ) 



3. DATA ANALYSIS 

3.1. H-a Spectroscopy 

We begin our analysis by determining the Ha equiva- 
lent width (EW) of all of our spectroscopic data, with no 
correction for the constant underlying photospheric ab- 
sorption line present. Following continuum normalization, 
the EWs of our Ritter data were computed using splot 
in IRAfI^. Equivalent width errors were computed using 
0-2 = N[hA/(S/N)]2(f,/f^), where N is the number of pixels 
across the Ha line, h^ is the dispersion in A^^, f, is the flux 
in the line, fc is the flux in the continuum, and S/N is the 
signal-to-noise ratio calculated in a line-free region of the 
spectra neighboring the Ha line ([Chalabaev fc MaillardI 
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Il983fl . For our HPOL data, the spectral regions 6350- 
6450 and 6700-6800A were used to constrain the contin- 
uum flux around Ha while the region 6520-6620A was used 
to measure the Ha line flux. Three nights of our tt Aqr 
HPOL dataset (MJD 2452860, 2452906, 2453181) sufl'ered 
from extremely poor wavelength calibration, requiring us 
to shift the three spectra from these nights by ■~80A. 

Our EW data are compiled in Table [1] (60 Cyg) and 
Table [2](7r Aqr). The EWs of our HPOL and Ritter ob- 
servations of 60 Cyg are systematically offset from one 
another. This is likely caused both by the different resolu- 
tion of the instruments used and by the continuum in our 
low-resolution HPOL spectra not being perfectly flt via 
the linear function derived from our 2 continuum fllters. 
We stress however that this constant offset has no effect 
on our analysis or conclusions. Comparison of four nights 
in which both instruments obtained spectroscopic obser- 
vations indicates this offset is ~1.9A; for plotting purposes 
only (Figure [T|), we have added a constant offset of 1.9A 
to our HPOL EWs measurement. Note that no offset has 
been added to the tabulated measurements in Table [T] 

3.1.1. 60 Cyg 

Figure [T] clearly demonstrates that 60 Cyg experi- 
enced a fundamental change at Ha throughout our 14 
years of spectroscopic monitoring data, during which the 
line changed fr o m ab sorption to emission to absorption. 
iKoubskv et al.l (|2000t ) qualitatively report a potentially 
similar event prior to our observations, but they did not 
publish their detailed EW measurements and only referred 
to the disk as being in "strong emission, an intermediate 
state, and pure absorption", thereby precluding a more 
comprehensive comparison of the events. Our data indi- 
cate that it took ~870 days for 60 Cyg's disk to transform 
from its strongest Ha emission state (1998 December 24; 
Table [1]) to its lowest state of pure absorption, the latter 
of which we interpret as the time when the disk had com- 
pletely dissipated (2001 May 9; Table[T|). Inspection of the 
decline state, seen in the right-panel of Figure [1] reveals 
that to flrst order the disk's emission strength followed a 
monotonic decrease void of any dramatic outbursts. After 
reaching its low state in 2001 May, 60 Cyg remained at a 
moderately low state characterized by occasional minor in- 
creases in the Ha EW through the last of our observations 
in 2006 December. 

Binary companions can influence the circumstellar 
disks of Be stars in several ways. They are be- 
lieved to_Jnflueiice_ the behavior of global disk oscilla- 
tions (lOktariani fc Ok azaki 2009) , truncate the outer disk 



(|Carciofi et al. 20091 ) . and potentially trigger mass-loss 



events from the centra l star for systems in very eccen- 
tric orbits (e.g. 5 Sco; iMiroshnichenko et al.ll200i[) . As- 
suming the ra nge of primary and se condary masses for 60 
Cyg given by IKoubskv et al] (|2000D and a circular orbit, 
60 Cyg's disk could be truncated at the tidal distance of 
^^120-135 Rs„„, assuming t he tidal radius is '--^0.9 R_Roc/ie 
(valid for small mass ratios: IWhitehurst fc King 199 ^ an d 
using the Roche radius approximation of lEggletonI (jl983l ). 
Aside from this possible truncation, we suggest it is un- 
likely that the binary companion is a fundamental con- 
tributor to 60 Cyg's disk- loss event. The time-scale of the 
disk-loss event (~870 days) corresponds to almost 6 com- 



plete orbits of the binary comp anion (assuming a 146.6 
day period; IKoubskv et al.|[2000l ). suggesting that the bi- 
nary does not influence the primary star (or its disk) in 
a manner similar to the h ighly eccentric 5 Sco system 
(IMiroshnichenko et al1l200H ). 

3.1.2. TT Aqr 

The time evolution of the Ha EW of tt Aqr (Figure [2]) 
also provides clear evidence of a disk- loss event, at a higher 
sampling frequency tha n available in previous publications 
(|Biorkman et al.ll2002l ). Our data indicate that it took at 
least ~2440 days for tt Aqr's disk to transform from its 
strongest state observed during our time coverage (1989 
October 26; EW = -24.6 A Table H]) to the hkely start of 
its minimum state of pure absorption (1996 July 4; EW = 
4.1 A Tabled]). While the overaU behavior of tt Aqr's de- 
cline phase shows a mostly monotonic decline in emission 
strength (left panel. Figure [2]), similar to that seen for 60 
Cyg, there is clear evidence of two significant interruptions 
to this trend. The decline of tt Aqr's Ha emission experi- 
enced a ^220 day flattening event beginning on 1991 May 
12 (JD 2448389) and a -200 day event beginning on 1994 
June 5 (JD 2449509), during which its EW not only flat- 
tened but exhibited evidence of temporarily increasing in 
emission strength (right panel. Figure H]). 

The time-scale for tt Aqr's disk-loss episode (—2440 
days) corresponds to roughly 29 complete orbits of its 
kn.own bi nary c ompanion (assuming a 84.1 day period; 
iBiorkman et al]|2002t ). which suggests that the passage of 
the binary was not responsible for triggering the disk-loss 
episode. Assuming the ra nge of primary and sec ondary 
masses for tt Aqr given by IBiorkman et al.l ()2002t ) and a 
circular orbit, we calculate that tt Aqr's disk could be trun- 
cated at the tidal distance of 80-89 l\.sun, assuming the 
tidal radius is —0.9 RRoche ()Whitehurst fc Kins'' 1 9|9lr) an d 
using the Roche radius approximation of EgglctonTi 19831 ) . 
Aside from this possible truncation, we suggest it is un- 
likely that the binary companion influences the primary 
star (or its disk) i n a manner similar to the h ighly eccen- 
tric 5 Sco system (jMiroshnichenko et al.ll200ll ). 

After the likely start of its minimum "pure absorption" 
state in 1996 July, tt Aqr's Ha EW remained generally 
unchanged until 2001 November, at which time a minor 
fiUing-in of the photospheric absorption line occurred and 
lasted until 2002 June. A similar strengthening event ap- 
peared to begin at or around 2003 October and persisted 
through our final observation on 2004 October, though 
the Ha emission strength never approached the values ob- 
served during the star's previous strong-disk phase. 

3.2. Total V-hand Polarization 

Theoretical models predict that the predominant phys- 
ical locations in the disk at which Ha photons are emit- 
ted and at which photons are scattered to produce the 
observed polarization differs in Be stars, with the polar- 
ization produced in the very inner disk region and Ha 
photons coming from throughout a much more extended 
radia l distance ([Carciofi fc Biorkmanll2066t ICarciofi et al.l 
I2OOI . For example, application of th ese models to ob- 
servat i onal data for the Be star C, Tau (jWisniewski et al.l 
I2007at IStefl et al.ll2009f) suggest that while the scattering 
events producing polarized photons are produced within 
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a few stellar radii, Ha photons are emitted over a sub- 
stantially broader region of the disk, extending out to ^40 
stellar radii. Therefore, our spectropolarimetric datasct 
provides unique simultaneous insight into the long-term 
behavior of the inner versus outer regions of 60 Cyg's and 
TT Aqr's circumstellar disks. To facilitate comparison with 
our Ha EW data, we have binned our spectropolarimetric 
data to replicate the Johnson V-band filter. 

3.2.1. 60 Cyg 

The time evolution of 60 Cyg's observed (intrinsic and 
ISP) V-band polarization is depicted in the top panels of 
Figure [T] and tabulated in Table [TJ While the observed 
magnitude of polarization mimics the similar long-term 
trend as the strength of the Ha line emission, the maxi- 
mum Ha emission lags the peak polarization observed on 
1998 August 24 by ~120 days, as depicted by the first 
two vertical lines in the right panels of Figure [TJ Simi- 
larly, the onset of the minimum Ha EW (third vertical 
line, Figure [1]) lagged the onset of the polarization low- 
state by '^190 days. Such a time lag between the on- 
set of polarimetric and spectroscopic features in Be stars 
has been previously reported i n the literature (see e.g. 
iPoeckert &: Marlboroughlll97 8a''b') and is indicative that 
the clearing of the disk proceeds in an "inside-out" man- 
ner. If we assume that the time period immediately fol- 
lowing the polarization maximum denotes the beginning 
of the disk-loss transformation, then the full time-scale for 
disk loss is simply the summation of the Ha EW decline 
time-scale plus the polarization time lag (~870 days -I- 
~120 days), ~1000 days. Thi s disk loss time-scale i s simi- 
lar to the ~700 day time-scale iPoeckert et al.l (jl982t ) noted 
for the loss of omicron And, using the same definition of 
the start and termination of the disk-loss phase. 

In Section 13.1.11 we noted the presence of minor, tem- 
porary increases in the Ha EW after 60 Cyg reached its 
diskless low state. Inspection of Figure [T] illustrates that 
one of these EW events near 2004 Aug 30-31 was accompa- 
nied by a ^0.3% V-band polarimetric outburst, strength- 
ening our interpretation that these minor outbursts repre- 
sent tenuous attempts of the system to reform a gaseous 
circumstellar disk. 

Though not plotted in Figure [1] for asthctic reasons, the 
PBO Lyot polarization observation from 1987 (Table [3]), 
which encompasses the V-band, is consistent with HPOL 
polarimetry obtained during epochs in which 60 Cyg was 
characterized by a strong disk. These da t a are also con- 
sistent with the results of iKoubskv et al] (|2000f) . who re- 
ported strong Ha emission present before (early 1980s) 
and after (later 1980s) our PBO Lyot observation. 

3.2.2. TT Aqr 

The time evolution of tt Aqr's observed V-band polariza- 
tion, depicted in the top panels of Figure [Hand tabulated 
in TableHl exhibits a substantially more complex morphol- 
ogy during the star's disk-loss episode as compared to 60 
Cyg. Archival V- band polarization pri or to the onset of 
our observations (jMcDavidI Il986l . Il99 0) , plotted as open 
circles in Figure [U indicate that the disk was in a stable 
high state between 1985-1987. Figure [2] clearly shows that 
the ~2440 day long decline of the Ha EW was interupted 
by two massive polarimetric outbursts. The ~390 day long 



first polarimetric event, demarked by the first and second 
vertical lines in the right panels of Figure [2l is clearly 
related to the '^220 day long first Ha EW event noted 
in Section [3.1.21 The ^730 day long second polarimetric 
event, demarked by the third and fourth vertical lines in 
the right panels of Figure ^ is similarily related to the 
^200 day second interuption of the Ha EW decline noted 
in Section 13.1.21 during which the EW briefly began to 
re-strengthen. 

Given the complex morphology of both the V-band po- 
larization outbursts and the interruptions of the Ha EW 
decline, it is more difficult to quantify the exact time lag 
between the onset of polarimetric features and the cor- 
responding effects in Ha. However, inspection of Figure 
[2] clearly shows that such a time lag does exist; for ex- 
ample, the total V-band polarization reaches its low state 
much earlier than the time at which the Ha EW reaches 
its low state. Similar to our interpretation of our 60 Cygni 
data, we deduce this time lag as evidence that tt Aqr's disk 
cleared in an "inside-out" fashion. 

In Section 13.1.21 we noted the presence of two minor, 
temporary increases in the Ha EW after the polarimetric 
data indicated that vr Aqr had reached its quiescent, disk- 
less state. We note that the final EW increase, which be- 
gan in 2003 October and was still ongoing at the termina- 
tion of our observational data in 2004 October, coincided 
with a corresponding temporary increase in the system's 
total V-band polarization, which peaked at 0.2% above 
the background interstellar polarization level on 2004 Au- 
gust 29. We observed no increase in the polarization cor- 
responding to the 2001 November - 2002 June Ha EW 
outburst, and speculate that we missed observing the on- 
set of this event due to the sparse time sampling of our 
polarimetric observations during this epoch. 

Though not plotted in Figure [T] for asthetic reasons, the 
PBO Lyot polarization observations from 1979-1980 (Ta- 
ble [J]), which encompass the V-band, are consistent with 
HPOL polarimetry obtained during epochs in which tt Aqr 
was characterized by a strong disk and also consistent with 
the general l evel of V-band polariza tion in the late 1970s 
presented bv lBiorkman et all (|2002D . 

3.3. Interstellar Polarization 

While polarimetry can be used as a unique diagnos- 
tic tool to investigate the properties of astrophysical ob- 
jects, proper interpretation of these data is inherently 
difficult, as observations can contain contributions which 
are both intrinsic and interstellar in origin. The identi- 
fication and complete removal of the interstellar polar- 
ization component (hereafter ISP), and a thorough un- 
derstanding of the accuracy of this removal, is thus a 
critical first step in interpreting pola rimetric data. As 
detailed by McLean fc Brownl (|1978D and discussed in 
lOuirrenbach et al. (199^, there are three techniques com- 
monly used to determine ISP along a line of sight: field 
star studies, emission line polarization analysis, and wave- 
length dependence studies. 

3.3.1. Field Star Technique 

Field star studies are the most basic of the ISP sep- 
aration techniques, and due to their large uncertainties, 
they are most useful when used in combination with other 
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methods. The technique involves averaging the observed 
polarization of stars at the similar distances to and small 
angular separations from the science target. It inherently 
assumes the existence of uniform interstellar medium con- 
ditions in the vicinity of the target and also assumes that 
the field stars themselves are devoid of any intrinsic po- 
larization component. While the technique can be highly 
successful when used to analyze the ISP properties of open 
clusters, where one can assume all members are located 
at the same distance and there are sufficient numbers 
of stars to lessen the effects of small intrinsic polariza- 
tion components in the field stars by b rute force averag- 
ing of a large number of stars (see e.g. IWisniewski et all 
l2007bf ). the apphcation of the techniq ue to isolated stars 
often yields imprecise ISP es tima tes (see McLean & Clarke 
119791 iBiorkman et aLlfToOSl and IWisniewski et al.ll2006l) . 

We were able to identify seven non-emission 
line stars with p o larim etric data in the catalog of 
iMathewson fc FordI (|1970[ ) spatially near the position of tt 
Aqr on the sky which appeared to be suitable field stars 
(see Table S]). Averaging these field stars yields the ISP 
parameters of Pism = 0.47% at 6ism = 130.1° . As we will 
show in Section I3.3.3[ this initial ISP estimate does not 
fully remove the interstellar polarization component from 
our data; hence, we do not adopt these values as our final 
ISP determination along the line of si ght t o tt A qr. 

Using the polarization catalog of HeilesI (|2000l ). we were 
only able to identify one non-emission line star with suit- 
able pol arimetric data haying a similar distance as 60 Cyg 
(418 pc, iPerrvman et al.l [19971 see Table g]). As seen in 
Figure [3l most of the stars spatially coincident with 60 
Cyg are in the foreground; moreover, the one star at a 
similar distance (250 - 550 pc) as 60 Cyg exhibits a signif- 
icantly different polarization. We conclude that the field 
star technique fails to yield a useful estimate of the ISP 
along the line of sight to 60 Cyg. 

3.3.2. H-a Line Depolarization 

If one assumes that the physical location of the emit- 
ting region of emission lines in an astrophysical ob- 
ject lies exterior to the region p roducing any intrin- 
sic po larization, as suggested by iHarrington fc Collins! 
(|1968[ ). one can measure the polarization across emis- 
sion lines such as Ha and use these values as measures 
of the ISP along the line of sight. While this tech- 
nique has bee n succ essfully employe d for a wide vari- 
ety of objects (iNord sieck et al. 200 11: iMever et al.l l2002t 
IWisniewski et al.l l2006f) . IMcLean fc Clarkd (|1979l ) note 
that the semi-empirical relationship used to isolate ISP 
via this technique is questionable and s ubject to nu- 
merous complications. iQuirrenbach et all (|1997( ) remark 
that the technique's assumption of a constant intrin- 
sic position angle across Ha is also doubtful and note 
that their own spectropolarimetri c observations s h ow Ha 
is not completely unp olariz ed. Poeckert et al.l (I1979D . 



lOudmaiier et all (|2005D . and IHarrington fc Kuhil l|2009t ) 
also demonstrate that Be stars can exhibit intrinsic line 
polarization features. The reliability of this technique is 
therefore uncertain at best. 

To apply this technique to our data, we first created 
artificial filters which sampled regions within the Ha line 
core (6520 -6620A) and continuum on both sides of the line 



(6350 -6450 A and 6700 -6800 A). We then determined the 
polarization within the Ha line and its nearby continuum, 
using these filters, for all data which showed clear spectro- 
scopic evidence of having Ha in emission. We plotted the 
continuum and line polarizations for each observation of 
TT Aqr and 60 Cyg on separate Stokes Q-U diagrams, and 
looked for the standard evidence of full line depolarization 
effects. Careful inspection of these figures revealed that 
the lines connecting each continuum-line observation set 
pointed in a different direction and that the Ha line polar- 
ization measurements were not clustered around one cen- 
tral location. These trends indicated that intrinsic polar- 
ization within the line was likely contaminating our data, 
thereby precluding our use of Ha line depolarization as an 
ISP constraint for both tt Aqr and 60 Cyg. 

3.3.3. Serkowski Law Fitting 

The wavelength dependence of interstellar polariza- 
tion can be exp ressed by the empirical Serko wski law 
(jSerkowski et al.l [1975. ) . as modified by Wil king et al.l 
([1981 : P(A) = P^ax exp K {\r,ax / \)\ ■ In the ab- 
sence of additional constraints, one can attempt to deter- 
mine the ISP along the line of sight to an object contain- 
ing both interstellar and intrinsic polarization if one has a 
priori knowledge of the positio n angle of the instrin sic po- 
larization component (see e.g. iPoeckert et al.|[T979f ). Our 
TT Aqr and 60 Cyg spectropolarimetric data are particu- 
larily unique in that they span an era in which each star 
enters a "normal B star" phase in which all observational 
evidence of the presence of a disk, including an intrinsic 
polarization component, had disappeared. The remnant 
polarization in these low-state data can be assumed to be 
fully interstellar in origin, thereby providing us the oppor- 
tunity to accurately fit their wavelength dependence using 
a modified Serkowski law. 

We begin by plotting our entire V-band datasets for 60 
Cyg and tt Aqr on a Stokes Q-U diagram as seen in Figures 
m and O repectively, where 



g = Pcos{29) 
U ^ Psin{2e) 



2\Q.5 



P= iQ' + U^) 
9 OMan-\U/Q) 



(1) 
(2) 

(3) 

(4) 



The time evolution of the polarimetric data in each of these 
figures follows a strict linear trend, indicating a constant 
scattering angle (i.e. position angle) regardless of the mag- 
nitude of polarization, as is expected for electron scatter- 
ing in a gaseous circumstellar disk. The slope of the best 
fit linear regression to our 60 Cyg data (Figure H]) was 
0.71 ±0.02. This indicates that the intrinsic polarization 
position angle for 60 Cyg is 9^, — 107.7 ±0.4°, and under 
the typical assumption that the maximum polarization oc- 
curs for 90° scattering it also implies that 60 Cyg's disk 
is oriented at a position angle of Odisk — 17.7° on the sky 
(measured in the convention of North to East). The slope 
of the best fit linear regression to our tt Aqr data (Figure 
O was -0.50 ±0.01. This translates into an intrinsic po- 
larization position angle, 6** of 166.7 ±0.1° for tt Aqr, and 
implies that tt Aqr's disk is oriented at a position angle of 
Odisk = 76.7° on the sky. 
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Using an iterative process, we next carefully examined 
the polarimetric data which corresponded to time-periods 
in which Ha was nearly or completely in absorption to 
help identify all observations which were devoid of an in- 
trinsic polarization component. Data comprised solely of 
an interstellar polarization component, which we labeled 
as "L" for polarization low-state in Tables [1] and [2l clearly 
correlated with one another in a clump on the Stokes Q- 
U diagram (see e.g. Q ~0% and U~0.1% in Figured]). 
Recall from Section 3.2 that H-a emission probes a much 
more extended region of a Be star's circumstellar environ- 
ment than optical polarization measurements; therefore, 
our definition of a polarization low-state only establishes 
that the density of material in the inner disk region dra- 
matically dropped (or was completly cleared) and there- 
fore does not demand that the outer circumstellar envi- 
ronment be similarly completely cleared of disk material. 
Since the intrinsic polarization arising from classical Be 
stars follows a tight linear trend, we can effectively isolate 
the presence of any residual intrinsic polarization in our 
polarization low-state data by rotating all of our obser- 
vations by the negative value of the intrinsic polarization 
position angle. 

For 60 Cyg, we rotated our data by -107.7°, thereby 
placing any small amount of residual intrinsic polarization 
in the Stokes Q parameter (where the prime denotes the 
rotated frame) and fully isolating the perpendicular com- 
ponent of the interstellar polarization in the Stokes U pa- 
rameter. After co-adding all low state polarization data in 
this rotated frame, we then fit a modified Serkowski law 
to the U data using the IDL subroutine mpcurvefit, as 
demonstrated in Figure [6l which plots the U component 
of the 60 Cyg data binned to a constant (Poisson photon 
statistic) polarization error of 0.01%. 60 Cyg's perpen- 
dicular ISP component (U data) is clearly well fit by a 
modified Serkowski law having parameters listed in Table 
\5\ There is very little signal in 60 Cyg's Q component 
(Figured]); hence, we could not fit the data with a modi- 
fied Serkowski law to extract the parallel ISP component. 
Instead, we overplotted our perpendicular ISP component 
on our low-state polarization data in the Stokes Q -U pa- 
rameter space and estimated the lower-limit of the parallel 
ISP component needed to ensure that our lowest state dat- 
apoint was purely interstellar in nature (see dashed line in 
Figure [7]). Our lower limit estimate of the parallel ISP is 
tabulated in Table along with the total ISP estimate 
derived from vectorally adding the perpendicular and par- 
allel ISP components. We find the total ISP for 60 Cyg 
to be characterized by Pmax = 0.112% ± 0.001%, Xmax — 
6977 ± 24lA, PA = 41.5°, and K = 1.17. The wavelength 
dependence of our polarization low state 60 Cyg data after 
subtraction of this total ISP, seen in Figure [5] binned to a 
constant polarization error of 0.015%, is centered about Q 
= 0% and U = 0% and exhibits no significant structure, 
thereby confirming the validity of our ISP determination. 

We employed a similar technique to determine the ISP 
for 77 Aqr. We first rotated our data by -166.7° to place 
any small level of residual intrinsic polarization in our low- 
state data in the Stokes Q parameter and isolating the 
perpendicular ISP component in the Stokes U parameter. 
We were able to obtain a good fit of a modified Serkowski 



law to the co-added U data again using mpcurvefit, as 
demonstrated in Figure [21 which plots tt Aqr's U compo- 
nent binned to a constant (Poisson photon statistic) po- 
larization error of 0.01%. We do note that a small fraction 
of the binned U data at the extrema of our wavelength 
coverage, e.g. below ^4000 A and above ~9500A, deviate 
from the best fit by <0.1%. These deviations are less than 
3-cr from the known instrumental systematic uncertain- 
ties (Section [5]); hence, we do not believe they represent 
real features. To confirm this, we also fit the data with 
an error weighted modified Serkowski-law, after adopting 
the known systematic instrumental uncertainties for the 
U- and I-band data, which were larger than the Poisson 
photon statistic errors. Similar to the result achieved for 
completely excluding the U- and I-band data from the fit, 
we found this exercise did not appreciably affect the resul- 
tant best fit ISP parameters. Since the number of binned 
data points which deviate from the fit is extremely small, 
relative to the total number of binned data points used the 
fit, these two modifications to the fitting routine had no 
effect on the resultant best fit. The Serkowski law param- 
eters of TT Aqr's perpendicular ISP component are listed 
in Table [i 

For the parallel ISP component, we included an addi- 
tional constant term in our modified Serkowski law fitting, 
"C" , to accomodate the possibility that a small level of in- 
trinsic polarization was present in our low state data and 
fixed the value of Xmax to that found in the perpendic- 
ular ISP component (Table [5]) , as this parameter should 
be the same for both the perpendicular and parallel ISP 
components. The resultant best fit to the Q component 
agrees reasonably well with the data, as seen in Figure [SI 
with small deviations present below ^^4000 A and above 
-9500A. These <0.05% deviations from the best fit are 
less than 3-ct from the known instrumental systematic un- 
certainties (Section [2]); hence, we do not believe they are 
real and, as for the perpendicular ISP fit, they do not sig- 
nificantly influence the best fit parallel ISP parameters. 
We find the total ISP for tt Aqr to be characterized by 
Pmax = 0.514% ± 0.001%, Xmax = 4959 ±15 A, PA = 
108.7°, and K — 0.83, as computed by vectorally adding 
the parallel and perpendicular ISP components listed in 
Table [5] The wavelength dependence of our polarization 
low state TT Aqr data after subtraction of this total ISP, 
seen in Figure [TU] binned to a constant polarization error 
of 0.005%, exhibits only a small amount of residual polar- 
ization (< 0.1%) in either the Stokes Q or U components 
and exhibits little structure outside of the aforementioned 
regions in which the HPOL CCD sensitivity is low. 

3.3.4. Comparison to Previous ISP Estimates 

While we are not aware of any successful attempts to 
char acterize the ISP alon g the line of sight to 60 Cyg (see 
e.g. iPoeckert et al.|[T979t ). two groups of authors have at- 
tempt ed to derive the ISP alo ng the line of sight to tt 
Aqr. [McLean fc Clarkd (|1979D estimated the ISP to be 
Pism ~ 0.45% and 9ism ~ 130° from analysis of field stars, 
Pism > 0.6% based upon tt Aqr's E(B-V) excess, and 
Pism < 0.4% and 9ism < 130° from studying the wave- 
length dependence of broad-band observations. The au- 
thors also note that the latter parameters were supported 
by Ha and H/3 line polarimetry. From these various meth- 
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ods. lMcLean k Clarkd (|l979t) adopted the final ISP values 
oiPmax ~ 0.36%, Xmax = 5400A, 0i ,m. = 120°, and placed 
no formal errors on these estimates. iPoeckert et al.l|i979t ) 
analyzed a heterogeneous mix of optical intermediate-band 
filter polarimetry, Ha polarimetry, and IR polarimetry to 
derive three different ISP estimates with large reduced 
values r anging from 5.6 - 10.7. While the Xmax value de- 
rived by IPoeckert et al.l (|1979[ ] agrees with the value we 
have derived (Table [S]) their PA and Pmax values differ by 
'^7° and up to 0.1% respectively. 

4. DISCUSSION 
4.1. Disk-Loss and Disk-Renewal Time-scales 

Identifying the mechanisms responsible for the disk-loss 
and disk-renewal phases of classical Be stars (i.e. the 
transition to and from a "Be" phase and a normal B-star 
phase) can be aided by an examination of the time-scale 
of some or all of the four main phases involved: the disk 
growth phase, the steady-state disk phase, the disk-loss 
phase, and the steady-state quiescent (diskless) phase. 

Our data poorly constrain the time-scale for the com- 
plete regeneration of an active disk state from a quiescent 
diskless state. Owing to the poor temporal sampling of 
our 60 Cyg data between 1992 and 1997 (Figure[T]), we can 
only place an upper limit on the regeneration of the active 
disk state observed in our data of <5 years. We do see 
evidence of temporary, correlated enhancements in both 
the V-band polarization and Ha EW during 60 Cyg's qui- 
escent, diskless state, including one moderate event which 
occurred between June and September 2004 in the V-band 
polarization data. Similarily, our long-term observations 
of TT Aqr do not cover a full disk regeneration cycle which 
would enable us to constrain this time-scale. However, 
we do observe one major, correlated enhancement in both 
the V-band polarization and Ha EW during the quies- 
cent phase (Figure 12) , again during the June - September 
2004 time-frame in the V-band polarization data. The 
sparse temporal sampling of 60 Cyg's and tt Aqr's out- 
burst episodes during their quiescent phases only enables 
us to place upper limits on the time-scales of their duration 
of < a few months. 

Ha monitoring of the B3 star fi Cen, revealed similar 
types of outburst events occuring during the star's qui- 
escent phase, which Hanuschik ct al. (1993) interpreted as 
evidence of individual blobs of material being injected into 
the star's inner circumstellar environment on time-scales 
of 2-5 days, and then viscously accreting back onto the star 
and/or decreting on time-scales of 20-80 days. Given the 
similarity in time-scales between fi Cen's quiescent-state 
Ha outburst events and the duration of the polarization 
component of quiescent state outbursts observed in 60 Cyg 
and TT Aqr, we propose that our data might be diagnos- 
ing one (or several) injection events which subsequently 
viscously dissipate, perhaps arising from the same type of 
non-radial pulsati on events thought to drive the injection 
events on /i Cen (jRivinius et al.lll998D . Such an interpre- 
tation is consistent with p revious suggestions of non-radial 
pulsations in bo th tt Aqr (jPeters fc Gie£|[2005l ) and 60 Cyg 
(|Koubskv et"ani2000D . 

As discussed in Sections 13.1.11 and 13.1.21 our data pro- 
vide excellent temporal coverage of the disk-loss events 
of both 60 Cyg and tt Aqr. 60 Cyg experienced a gen- 



erally monotonic reduction in Ha emission strength from 
its maximum to minimum state over a time-scale of ~870 
days. We interpret the '^120 day lag of the maximum in 
the Ha EW from the maximum V-band polarization as 
evidence that the disk cleared in an "inside-out" man- 
ner. The total dissipation time-scale of 60 Cyg's disk, 
e.g. from the maximum polarization level to the mini- 
mum Ha EW level, of ~1000 days is similar to that ob- 
served fo£_the_dissi£ationo£^^ And's disk (~700 
days; IClark et al.ll2003f) . IClark et all ([2003) noted that 
omicron And's disk dissipation was consistent with the 
expected viscous tim e-scale of the disk, assumin g the vis- 
cosity parameter a (jShakura &: Svunvaevlll973[) is ~0.1. 
Assuming that 60 Cyg's disk also dissipated on a viscous 
time-scale, we can derive the viscosi t y par ameter a using 
equation 19 of iBiorkman fc Carciofil |2005( ): 

td^ffus^on = (0.2yr/a) * (r/i?,)°-^ (5) 
. Assuming that the bulk of the scattering events produc- 
ing the o bserved polarization occurs at a radial distance of 
^2.5 R* (jCarciofi et al.|[2009,) and a ^-^800 day time-scale 
for the polarization to transform from maximum to mini- 
mum state (see Figure [1]), we derive a viscosity parameter 
of ^^0.14 for 60 Cyg, which is similar to that found for 
omicron And. By contrast, the loss of tt Aqr's disk as 
measured by its Ha EW proceeded a factor of at least 2.4 
times slower (e.g. ~2440 days) than 60 Cyg's disk loss, 
and Figure [2] clearly illustrates the presence of two ex- 
tended ('^390 day and '^730 day) outburst events which 
temporarily stalled the long-term disk loss event. Like 
60 Cyg, TT Aqr also exhibited a lag between onset of the 
low-state of the V-band polarization and the low-state of 
the Ha EW, indicating the disk cleared in an "inside-out" 
manner. The exquisite temporal sampling of our tt Aqr 
data suggests that temporary increases in the mass injec- 
tion rate into the inner disk may produce substantial de- 
viations from a monotonic, viscous dissipation of Be disks. 
We therefore suggest that claims of inconsistencies in the 
formation versus dissipation time-scales of disks with a 
standard viscous disk model ba s ed on sparsely sampled 
data, as made bv IVinicius et al.l ()2006[ ). might be caused 
by under-sampling the type of short-lived outbursts ob- 
served in our tt Aqr data. 

4.2. Temporal Changes in the Intrinsic Polarization 
Position Angle 

While our total V-band polarization data for both 60 
Cyg and tt Aqr both clearly followed a linear trend when 
plotted on a Stokes Q-U diagram (Figures S] and O, close 
inspection of these figures reveals tentative evidence of 
small-scale deviations in some of the data from each star's 
intrinsic polarization position angle, which seems to in- 
crease with the strength of the intrinsic polarization com- 
ponent. To quantify this effect, we calculated the absolute 
value of the error weighted deviance of each data point 
from the best fit line denoting the star's intrinsic PA. We 
plot these values as a function of the magnitude of the in- 
trinsic polarization component for both 60 Cyg (top panel) 
and TT Aqr (bottom panel) in Figure lll[ where filled tri- 
angles represent polarization low-state data (Tables [1] and 
(21) and open diamonds represent high-state data. 

Several trends are immediately apparent. For both 60 
Cyg and tt Aqr, the bulk of the data exhibit a < 6tT residual 
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from the bulk value of the intrinsic PA value; moreover, 
very few data points at low intrinsic polarization levels 
deviate at the > 6a level. Just 2 of 15 (13%) of 60 Cyg 
data below an intrinsic polarization level of 0.3% exhibit 
a > 6(T deviation while a mere 4 of 59 (7%) of tt Aqr 
data below an intrinsic polarization level of 0.6% exhibit 
a > 6cr deviation. By contrast, a substantial, statistically 
significant number of data at larger intrinsic polarization 
levels exhibit large deviations from the bulk intrinsic PA. 
8 of 15 (53%) of 60 Cyg data above an intrinsic polariza- 
tion level of 0.3% exhibit a > 6a deviance while 15 of 32 
(47%) of TT Aqr data above an intrinsic polarization level 
of 0.6% exhibit a. > 6a deviance. If these trends are in fact 
real, they would indicate that disk material has a stronger 
likelihood to be located slightly above or below the disk 
midplane when the inner disk density is larger (i.e. when 
the intrinsic polarization component is larger). For exam- 
ple, the position angle of the tt Aqr data point exhibiting a 
deviance of ^26a deviated by 2.43° from the bulk intrinsic 
PA of the disk of of 166.7°. 

Before continuing our interpretation of this phe- 
nomenon, we consider the appropriateness of our choice 
of a deviance cutoff, above which we believe the data ex- 
hibit a bona-fide deviation from the intrinsic polarization 
PA (e.g. 6a in Figure [TT|) . Merely setting a threshold 
deviance value above which everything is considered real 
could lead to several false discoveries, owing to random 
variations in the data. The false- d iscove ry rate (FDR) 
me thod, develope d by Miller et"al] (l200lh ari d emp loyed 
by iHopkins et all (|2002[ ) and iKowalski et"al] (|2009[ ). al- 
lows us a means to quantify the percentage of false posi- 
tives above our choice of a threshold deviance value. We 
model the null distribution (random noise) using the de- 
viance values below an intrinsic polarization level of 0.6% 
for TT Aqr and 0.3 % for 60 Cyg. Employing the IDL routine 
in Appendix B of lMiller et al.l (j2001f) . we determined that 
setting a threshold to the deviance value of 7.687 for tt Aqr 
predicts a false-discovery rate of 32%. In other words, of 
the 19 discoveries above a deviance of 7.687 (spanning the 
entire intrinsic polarization range in Figure [TT]), no more 
than 32% are caused by spurious deviations. If we assume 
that the four discoveries below an intrinsic polarization of 
0.6% are in fact false discoveries, this implies that only 
2 of the discoveries above 0.6% are false, e.g. 13 of the 
15 data points above the 6a level in the bottom panel of 
Figure [TT] are real. Even in the worst case scenario that 
the entire 32% of the sample of false discoveries resides 
in the discoveries we claimed above the 0.6% intrinsic po- 
larization level, a full 68% of this sample would still be 
real, statistically sound detections of deviations from the 
bulk polarization position angle. Similarily, for 60 Cyg, a 
false-discovery rate of < 33% sets the threshold deviance 
value at > 6.52; hence, of the 9 observed values above this 
deviance level, at most 3 could be false discoveries. Based 
on this FDR analysis, we conclude most of our claimed 
deviations from the bulk PA of both 60 Cyg and tt Aqr 
are in fact real. 

To further explore the origin of these PA deviations, we 
examine the deviance values computed for Figure [11] as a 
function of time and compare the trends in these data to 
the temporal evolution of the V-band polarization. Sev- 
eral clear trends are seen in the data for tt Aqr (Figure [T^ : 



1) a substantial number of deviance values exceeding the 
6a level coincided with the latter-portion of the ^390 day 
polarimetric outburst noted in Section [3.2.21 and 2) there 
was no corresponding increase in the PA deviance coin- 
ciding with the ~730 day polarimetric outburst reported 
in Section 13.2.21 Inspection of the right panel of Figure 
[T2l which depicts the behavior of the PA deviance during 
the ^390 day outburst, reveals tentative evidence of two 
separate components in the deviance. The rise of the first 
component of the deviance outburst (points A to B; Fig- 
ure [T^ spanned < 54 days while the decline of this first 
component (points B to C; Figure [T^ spanned < 13 days. 
The time-scale of the second component exhibited similar- 
ities, with a rise time (points D to E; Figure [T^ of < 35 
days and a decline time (points E to F; Figure [T^ of < 10 
days. We note that the frequency of our observations do 
not trace the decline phase at a sufficient sampling rate to 
diagnose the precise time-scale of the decline; rather, the 
time-scales we quote are merely upper limits. We can not 
exclude the possibility that the PA deviations returned to 
normal faster than the 10-13 days we quote. 

A similar trend is observed in the temporal evolution 
of the PA deviance of the 60 Cyg data (see Figure [T3]) . 
A large outburst is seen in the deviance which coincides 
with the epoch of the single polarimetric outburst noted 
in Section 13.2.11 Although the sampling of our 60 Cyg 
polarimetry is much more coarse compared to tt Aqr, Fig- 
ure [T2] illustrates that the < 26-day rise of the deviance 
outburst (points A to B) and subsequent < 14-day decline 
(points B to C) qualitatively matches the rise and decline 
time-scales observed for tt Aqr. 

4.2.1. PA Changes: Evidence of Structure in the Inner 
Disk or Blob Injections? 

Variability in the intrinsic polarization position an- 
gles of classical Be stars is not typically expected for 
a sta ndard axisymmetric disk system, al though short- 
term ( Carciofi et al.l 2007f) and long-term (jCarciofi et al.l 
2007t T^naka et al.l 12007 ) changes have been reported. 
Carciofi et al.l ( 20071 ) attributed variability on time-scales 
of hours to weeks in the Achernar system as evidence of the 
injection of blobs into the inner disk, leading to changes 
in the observed level of polarization and small changes in 
the polarization position angle which persist on order sev- 
eral wee ks until the bl o b is c ircularized into a ring. By 
contrast, iTanaka et al.l (|2007l ) attribute the decade time- 
scale evolution of Pleione's intrinsic polarization PA to 
disk precession caused by the system's binary companion. 
We interpret the short-term variability (tens of days) of 
the intrinsic polarization position angle of both tt Aqr and 
60 Cyg as evidence of low-amplitude, temporary depar- 
tures from a simple axisymmetric disk density within the 
inner disk regions (several stellar radii) of these systems. 

The deviations from axisymmetry we observe could be 
diagnostics of the injection of new blobs into the inner disk 
region, which subsequently circularize into a ring, thereby 
acting as an effective w a rp of the inner disk (see e.g. Fig- 
ure [14]). ICarciofi et al.l (|2007t ) reported variability in the 
polarization magnitude and position angle of the Be star 
Archernar, and offered a similar interpretation. These au- 
thors noted that the circularization of a single blob would 
occur on time-scales of days and that the time-scale for 
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a ring to incorporate into the inner disk region would be 
several weeks. The newly injected material would then be 
re-accreted onto the central star if the mass-loss episode 
was a singular event, whereas a steady stream of mass-loss 
episodes w ould cause the disk size and density to grow 
over time (jCarciofi et alj [20071 ) . As the PA variations in 
the events observed in both tt Aqr and 60 Cyg grew over 
time-scales of 1-2 months and decayed in < 1-2 weeks, we 
speculate that perhaps the departure from axi-symmetry 
we witnessed was caused by the injection of multiple blobs 
into these systems which circularized into a ring and were 
later incorporated into the inner disk region. 

Alternatively, the observed deviations from axisymme- 
try could be caused by a blob or series of blobs launched 
from a non-equatorial latitude from the stellar photo- 
sphere (Figure [Ti)) . Such blob(s) would still circularize on 
time-scales of days to weeks, but owing to the latitude of 
their ejection they would circularize into an inclined (non- 
coplanar) orbit with respect to the sem-major axis of the 
bulk of the disk material. One could imagine this inclined 
ring eventually incorporating itself into the plane of the 
bulk of the disk material owing to viscous drag (or other) 
forces. Thus we speculate that perhaps our observation 
of a temporary deviation in the polarization position an- 
gle might be diagnosing the presence of a newly injected 
blob, which is subsequently circularized, at an inclined or- 
bit with respect to the plane of the pre-existing disk ma- 
terial. We plan to explore the feasibility of both of the 
above scenarios usin g models similar to those presented in 
ICarciofi et aP (I2007D in a future publication. 

We note that one of the polarimetric outbursts of /pi 
Aqr, occuring in 1994 June, was not accompanied by a 
corresponding variation in the polarization position an- 
gle (see e.g. Figure [T^ . Unlike the other outbursts, our 
temporal coverage of the 1994 event was especially sparse. 
TT Aqr's V-band polarization was elevated in observations 
obtained on 5 June 1994, 7 June 1994, and 22 June 1994. 
However, the onset of this outburst is unkown, as no obser- 
vations were obtained between 10 September 1993 and the 
elevated 5 June 1994 observation. Similary, there was a ^6 
week gap between the last elevated observation on 22 June 
1994 and the subsequent observation on 3 August 1994, 
which merely places an upper limit on the time-scale of the 
polarization decline (of < 6 weeks). Given this sparse data 
sampling, we speculate that it is possible that we simply 
"missed" observing the injection events which would have 
produced variations in the PA, and instead only sampled 
the aftermath of these events in which newly injected ma- 
terial had already circularized in the disk. Alternatively, 
it is also possible that not all mass injection events yield 
deviations from axi-symmetry. Clearly, additional obser- 
vations of more polarimetric outburst events, at a higher 
sampling frequency, are needed to better constrain this 
phenomenon. 

4.3. Long-term Disk Loss/Disk Renewal Behavior of 60 

Cyg 

While our long-term Lyot and HPOL polarimetric mon- 
itoring of TT Aqr only covered 1 complete cycle of a Be to 
normal-B transition, our polarimetric observations of 60 
Cyg actually span multiple disk loss/disk renewal transi- 
tions. Specifically, there is clear evidence from the long- 



term behavior of both the V-band polariza tion and H-a 
EW (see Figure [2 also lKoubskv et al.|[200l that the cir- 
cumstellar disk diagnosed by our 1987 PBO Lyot polari- 
metric observation had dissipated by the early 1990s and 
that our polarimetric data from 1997-2000 was probing a 
physically new circumstellar disk. Interestingly, after sub- 
tracting the ISP from both the PBO Lyot and PBO HPOL 
polarimetry, we find the exact same (intrinsic) polariza- 
tion position angle in both data sets. This indicates that 
the older (late 1980s) and most recently generated (1997- 
2000) disks shared the same equatorial plane on the sky, 
and suggests a similar formation mechanism was respon- 
sible for both events. We note that this scenario is funda- 
mentally different than the misaligned double-d isk struc- 
ture reported for Pleione bv lTanaka et al.l ()2007D . whereby 
the authors suggested that the PA of the disk was precess- 
ing due to interactions with the star's binary companion. 

5. FUTURE WORK 

We have shown that the relatively rare occurence of a 
disk-loss episode in classical Be stars provides a unique 
opportunity to constrain the interstellar polarization com- 
ponent along the line of sight to these stars in exquisite 
detail. While we have analyzed the behavior of the total 
polarization of our dataset, our accurate ISP determina- 
tion will facilitate the improved removal of the interstel- 
lar component and thereby enable the time evolution of 
the intrinsic polarization component of these data to be 
studied in detail. We are currently analyzing the behav- 
ior of the intrinsic polarization of these stars and plan to 
present results of detailed modeling of these syste i ns, us - 
ing the Monte Carlo code of ICarciofi fc BiorkmanI (j2006[ ). 
in a future publication. 

We also encourage observers to follow-up on the results 
presented in Section 14.21 which indicate that the polar- 
ization PA of classical Be star disks may exhibit small- 
scale deviations during polarization outbursts. We have 
suggested that the origin of these deviations might be ev- 
idence of newly injected blobs in the disk which have not 
yet circularized into a ring, or evidence of new blobs which 
were launched from a non-equatorial latitude thereby pro- 
jecting them into an inclined circular orbit with respect to 
the plane of the pre-existing disk matrial. We encourage 
high time cadence observations which diagnose the inner 
regions of such disk systems, such as polarimetry and/or 
infrared spectroscopy, to further constrain the origin of 
this phenomenon. 

6. SUMMARY 

We have presented the results of 127 spectropolarimetric 
observations of the Be star tt Aqr spanning 15 years, and 
35 spectropolarimetric and 65 Ha spectroscopic observa- 
tions of the Be star 60 Cyg spanning 14 years. These data 
trace the evolution of each star as it transitions from hav- 
ing a gaseous circumstellar disk ( "Be phase" ) to a diskless 
state ("normal B-star phase"). We find: 

• 60 Cyg's disk emission strength, as measured by its 
Ha EW, declined from maximum emission to pure 
absorption monotonically in ~870 days. The max- 
imum observed Ha EW lagged the maximum ob- 
served V-band polarization by ~120 days, indicating 
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that the total disk dissipation time-scale was ~1000 
days (~870 + ~120), and that the disk clearing pro- 
ceeded in an "inside-out" fashion. We found the dis- 
sipation of disk, as measured from the time-scale of 
the decrease in the observed polarization, occured 
over a viscous time-scale if we adopted a viscosity 
parameter a of 0.14. 

• TT Aqr's disk emission strength as measured by its 
Ha EW declined from maximum emission to pure 
absorption over a time-frame of ~2440 days. Two 
extended outburst events lasting ^390 and ~730 
days respectively were observed during tt Aqr's disk- 
loss phase, which had the effect of temporarily 
stalling the loss of the disk. An observed lag in 
the onset of the minimum Ha EW strength from 
the minimum V-band polarization indicates that tt 
Aqr's disk clearing also proceeded in an "inside-out" 
fashion. 

• The time-scale of the disk-loss events in 60 Cyg and 
TT Aqr corresponds to almost 6 and 29 complete or- 
bits of each star's binary companion respectively. 
This suggests that each star's binary companion 
does not influence the primary star (or its disk) in a 
manner similar to the highly eccentric S Sco system 
(Miroshnichenko et al. 2001). 

• We fit the wavelength dependence of multi-epoch 
spectropolarimetric observations obtained during 
the quiescent diskless phase of each star with a mod- 
ified Serkowski law to derive the interstellar polar- 
ization along the line of sight to each star. We deter- 
mine the best fit modified Serkowski law parameters 
for 60 Cyg to be P„a. = 0.112% ± 0.001%, X„,ax = 
6977A, PA = 41.5°, and K = 1.17, while for tt Aqr 
we derive the parameters Pmax — 0.514% ± 0.001%, 
A„a^ = 4959A, PA = 108.7°, and K = 0.83. 

• We find that the position angle of intrinsic polariza- 
tion arising from 60 Cyg's disk is 0, = 107.7±0.4°, 
indicating that the disk is oriented on the sky at 
a position angle of Odisk — 17.7° (measured North 
to East). We find that the position angle of intrin- 
sic polarization arising from tt Aqr's disk is 6^ = 
166.7±0.1°, indicating that the disk is oriented on 
the sky at a position angle of Odisk = 76.7°. 

• We detect clear evidence of small outburst events 
in the V-band polarization and Ha EW during 
the quiescent diskless phase of both tt Aqr and 60 



Cyg, which persist for several months in the V- 
band polarization. As these events have similar 
time-scales as the blob injection events likely in- 
duced b^;_j]£n;Tadml_pulsations in the Be star j j, 
Cen (|Hanuschik et all Il993l : iRivinius et all Il998l ). 
we suggest these outburst events indicate the injec- 
tion of one (or several) blobs of material from the 
stellar photosphere which are subsequently accreted 
back onto the central star and/or decreted on a vis- 
cous time-scale. 

• We also detect evidence of deviations from the in- 
trinsic polarization position angle in both 60 Cyg 
and TT Aqr which persist for 1-2 months and coin- 
cide with large outbursts in the total V-band po- 
larization, and confirm these events are real via an 
examination of the false discovery rate for each ob- 
ject. For TT Aqr, we determine the magnitude of one 
of these deviations to be ~2.43° from the bulk in- 
trinsic polarization PA of the disk (166.7°). These 
results are indicative of deviations from a simple 
axisymmetric disk density structure in each system. 
We propose that either we are witnessing the injec- 
tion and subsequent circularization of new blobs into 
the inner disk reg i on in a similar manner as noted 
bv ICarciofi et al.l (|2007l ) or the injection and sub- 
sequent circularization of new blobs at an inclined 
orbit to the plane of pre-existing disk material (see 
Figure [Ti)) . 



We thank our referee, David Harrington, for providing 
helpful suggestions which improved the quality of this pa- 
per. We also thank Brian Babler for his invaluable assis- 
tance with various aspects of PBO HPOL data, the PBO 
and Ritter science teams for providing observing and data 
reduction support, and Kenneth H. Nordsieck for provid- 
ing access to the PBO Lyot and HPOL spectropolarime- 
ters. JPW thanks A. Carciofi and A Okazaki for pro- 
viding helpful comments regarding the effects of binary 
companions in Be star systems. JPW acknowledges sup- 
port from NSF Astronomy & Astrophysics Postdoctoral 
Fellowship AST 08-02230, and ZHD acknowledges partial 
support from the University of Washington Pre-MAP pro- 
gram. HPOL observations were supported under NASA 
contract NAS5-26777 with the University of Wisconsin- 
Madison. Observations at the Ritter Observatory are sup- 
ported by the NSF under the PREST grant AST 04-40784. 
This study has made use of the SIMBAD database, oper- 
ated at CDS, Strasbourg, France, and the NASA ADS 
service. 



REFERENCES 



Bjorkman, J. E., & Cassinelli, J. P. 1993, ApJ, 409, 429 
Bjorkman, K.S. et al. 1998, ApJ, 509, 904 

Bjorkman, K.S., Miroshnichenko, A.S., McDavid, D., & Pogrosheva, 

T.M. 2002, ApJ, 573, 812 
Bjorkman, J.E. & Carciofi, A.C. 2005, in ASP Conf Ser. 337, The 

Nature and Evolution of Disks Around Hot Stars, ed R. Ignace & 

K. Gayley (San Francisco: ASP), 75 
Brown, J.C., CassineUi, J.C., & Maheswaran, M. 2008, ApJ, 688, 

1320 

Carciofi, A.C. & Bjorkman, J.E. 2006, ApJ, 639, 1081 
Carciofi, A.C, Magalhaes, A.M., Leister, N.V., Bjorkman, J.E., & 
Levenhagen, R.S. 2007, ApJL, 671, 49 



Carciofi, A.C, Okazaki, A.T., le Bouquin, J-B., Stefi, S., Rivinius, 
Th., Baade, D., Bjorkman, J., & Hummel, CA. 2009, A&A, 504, 
915 

Cassinelli, J. P., Brown, J.C, Maheswaran, M., Miller. N.A., & Teller, 

D.C 2002, ApJ, 578, 951 
Chalabaev, A. & Maillard, J. P. 1983, A&A, 127, 279 
Clark, J.S., Tarasov, A.E., & Panko, E.A. 2003, A&A, 403, 239 
Clarke, D. & Bjorkman, K.S. 1998, A&A, 331, 1059 
Coyne, G.V. & Kruszewski, A. 1969, AJ, 74, 528 
Cranmer, S.R. 2005, ApJ, 634, 585 
Cranmer, S.R. 2009, ApJ, 701, 396 



11 



Doazan, V., Franco, M., Rusconi, L., Sedmak, G., & Stalio, R. 1983, 

A&A, 128, 171 
Eggleton, P.P. 1983, ApJ, 268, 368 

Hanuschik, R.W., Dachs, J., Baudzus, M., & Thimm, G. 1993, A&A, 
274, 356 

Harries, T.J., Bablcr, B.L., & Fox, G.K. 2000, A&A, 361, 273 
Harrington, J. P. & Collins, G.W. 1968, ApJ, 151, 1051 
Harrington, D.M. & Kuhn, J.R. 2009, ApJS, 180, 138 
Hcilcs, C. 2000, AJ, 119, 923 

Hopkins, A.M., Miller, C.J., Connolly, A.J., Genovese, C, Nichol, 

R.C., & Wasscrman, L. 2002, A,T. 123, 1086 
Koubsky, P. et al 2000, A&A, 356, 913 

Kowalski, A.F., Hawley, S.L., Hilton, E.J., Becker, A.C., West, A.A., 

Bochanski, J.J., & Scsar, B. 2009, AJ, 138, 633 
Lec, U., Osaki, Y., & Saio, H. 1991, MNRAS, 250, 432 
Lupie, O.L. & Nordsieck, K.H. 1987, AJ, 93, 214 
Mathewson, D. S., & Ford, V. L. 1970, Mem.R.A.S., 74, 139 
McDavid, D. 1986, PASP, 98, 572 
McDavid, D. 1990, PASP, 102, 773 
McDavid, D. 1994, PASP, 106, 949 
McDavid, D. 1999, PASP, 111, 494 
McLean, I. S., & Brown, J. C. 1978, A&A, 69, 291 
McLean, L S., & Clarke, D. 1979, MNRAS, 186, 245 
McSwain, M.V., Huang, W., & Gics, D.R. 2009, ApJ, 700, 1216 
Meyer, J.M., Nordsieck, K.H., & Hoffman, J.L. 2002, AJ, 123, 1639 
Miller, C.J. et al. 2001, AJ, 122, 3492 
Miroslmichenko, A.S. ct al. 2001, A&A, 377, 485 
Morrison, N.D., Knauth, D.C., MuUiss, C.L., & Lee, W. 1997, PASP, 

109, 676 
Neiner, C. 2002, A&A, 388, 899 

Nook, M. 1990, Ph.D. dissertation. The University of Wisconsin- 
Madison 

Nook, M., Cardelh, J.A., & Nordsieck, K.H. 1990, AJ, 100, 2004 
Nordsieck, K.H. et al. 2001, in ASP Conf Ser. 233, P Cygni 2000: 400 

Years of Progress, ed. M. de Groot & C. Sterken (San Francisco: 

ASP), 261 

Nordsieck, K. H., & Harris, W. 1996, ASP Conf. Ser. 97, Polarimetry 
of the Interstellar Medium, ed. W. G. Roberge & D. C. B. Whittet 
(San Francisco: ASP), 100 

Okazaki, A.T. 1991, PASJ, 43, 75 

Oudmaijcr, R.D., Drew, J.E., & Vink, J.S. 2005, MNRAS, 364, 725 
Oktariani, F. & Okazaki, A.T. 2009, PASJ, 61, 57 
Ferryman, M.A.C. et al 1997, A&A, 323, 49 



Peters, G.J. & Gies, D.R. 2005, in ASP Conf Ser. 337, The Nature 
and Evolution of Disks Around Hot Stars, ed R. Ignace & K. 
Gayley (San Francisco: ASP), 294 
Poeckert, R., & Marlborough, J. M. 1978a, ApJ, 220, 940 
Poeckert, R. & Marlborough, J.M. 1978b, ApJS, 38, 229 
Poeckert, R., Bastien, P., & Landstreet, J. D. 1979, AJ, 84, 812 
Poeckert, R., Gulliver, A.F., & Marlborough, J.M. 1982, PASP, 94, 
87 

Porter, J.M. & Rivinius, T. 2003, PASP, 115, 1153 

Quirrenbach, A., Bjorkman, K. S., Bjorkman, J. E., Hummel, C. A., 

Buscher, D. F., Armstrong, J. T., Mozurkewich, D., Ellas, N. M., 

& Babler, B. L. 1997, ApJ, 479, 477 
Rivinius, Th., Baade, D., Stcfl, S., Stahl, O., Wolf, B., & Kaufer, A. 

1998, A&A, 333, 125 
Scrkowski, K., Mathewson, D. S., & Ford, V. L. 1975, ApJ, 196, 261 
Shakura, N.I. & Syunyacv, R.A. 1973, A&A, 24, 337 
Slcttebak, A. 1982, ApJS, 50, 55 
Stefl, S. et al. 2009, A&A, 504, 929 
Struve, O. 1931, ApJ, 73, 94 

Tanaka, K., Sadakane, K., Narusawa, S., Naito, H., Kambe, E., 

Katahira, J., & Hirata, R. 2007, PASJ, 59, L35 
Townsend, R.H.D., Owocki, S.P., & Howarth, I.D. 2004, MNRAS, 

350, 189 

Undcrhill, A. & Doazan, V. eds. 1982, B Stars with and without 

Emission Lines (NASA SP-456; Washington, DC: NASA) 
Vinicius, M.M.F., Zorcc, J., Leister, N.V., & Lcvcnhagen, R.S. 2006, 

A&A, 446, 643 
Whitehurst, R. & King, A. 1991, MNRAS, 249, 25 
Whitney, B.A. & Clayton, G.C. 1989, AJ, 98, 297 
Wilking, B. A., Lebofsky, M. J., & Rieke, G. H. 1982, AJ, 87, 695 
Wisniewski, J. P., Babler, B.L., Bjorkman, K.S., Kurchakov, A.V., 

Meade, M.R., & Miroshnichenko, A.S. 2006, PASP, 118, 820 
Wisniewski, J. P., Kowalski, A.F., 13jorkman, K.S., Bjorkman, J.E., 

& Carciofi, A.C. 2007a, Ap.IL, 656, 21 
Wisniewski, J. P., Bjorkman, K.S., Magalhacs, A.M., Bjorkman, J.E., 

Meade, M.R., & Pcreyra, A. 2007b, ApJ, 671. 2040' 
Wolff, M. J., Nordsieck, K. H., & Nook, M. A. 1996, AJ, 111, 856 
Wood, K., Bjorkman, J. E., Whitney, B. A., & Code, A. D. 1996a, 

ApJ, 461, 828 

Wood, K., Bjorkman, J. E., Whitney, B. A., & Code, A. D. 1996b, 
ApJ, 461, 847 

Wood, K., Bjorkman, K. S., & Bjorkman, J. E. 1997, ApJ, 477, 926 



12 



Table 1 
Summary of 60 Cyg Data 



Date 


JD 


Obs 


Ha EW 


V-band Pol 


PA 


State 








A 


% 


degrees 




1992 Aug 3 


2448837 


PBO 


4.5 ± 0.1 


0.074 ± 0.006 


58.2 


H 


1992 Aug 31 


2448865 


PBO 


3.5 ± 0.1 


0.108 ± 0.004 


52.2 


H 


1997 Nov 10 


2450732 


PBO 


-5.5 ± 0.1 


0.631 ± 0.003 


102.1 


H 



Note. — A sample summary of spectroscopic and spectropolarimetric observations of 60 Cygni; the full version of this Table is available in 
the online-edition of this Journal. The third column denotes the location at which the data were obtained (PBO = Pine Bluff Observatory, 
RIT = Ritter Observatory), while the seventh column denotes a qualitative assessment of whether 60 Cyg is in a polarimetrically high state 
(H), indicative of the presence of a strong disk, or in a polarimetrically low state (L), in which evidence of the presence of a disk is absent 
(see Section l3.3.3l for further discussion of these designations). Note that we have labeled the two observations in 1992 August as being "high 
state" , primarily because they are isolated from subsequent observations by ~5 years. 
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Table 2 
Summary of tt Aqr Data 



Date 


JD 


Obs 


Ha EW 

A 


V-band Pol 

% 


PA 
degrees 


State 


1985 May 15 


2446200 


McD^ 




1.53±0.03 


158.8 




1986 Jun 15 


2446596 


McD^ 




1.47±0.02 


159.1 




1987 Aug 15 


2447022 


McD2 




1.50±0.03 


157.4 




1988 Jul 15 


2447357 


McD2 




1.04±0.06 


155.0 




1989 Jun 15 


2447692 


McD2 




0.97±0.04 


152.1 




1989 Aug 8 


2447747 


PBO 


-23.5±0.1 


1.277±0.005 


156.5 


H 


1989 Oct 26 


2447826 


PBO 


-24.6±0.1 


0.879±0.003 


148.9 


H 



Note. — A sample summary of spectroscopic and spectropolarimetric observations of vr Aquarii; the full version of this Table is available in 
the online-edition of this Journal. The third column denotes the location at which the data were obtained (PBO = Pine Bluff Observatory, 
McD = McDonald Observatory), while the seventh column denotes a qualitative assessment of whether tt Aqr is in a polarimetrically high 
state (H), indicative of the presence of a strong disk, or in a polarimetrically lo w state (L) , in w hich evidence of t he presenc e of a d isk is absent 
(see Section 13.3.31 for further discussion of these designations). ^ data are from lMcDawdl HTsM) . ^ data are from lMcDalodl digOOh . 
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Table 3 
PBO Lyot Polarimetry 



Target Star 


JD 


Date 


%Pol. 


P.A. 


Wavelength Coverage 








% 


degrees 


A 


60 Cyg 


2446998 


1987 Jul 22 


0.44 ± 0.02 


104 


4600-7200 


TT Aqr 


2444116 


1979 Aug 31 


1.17 ± 0.03 


153 


4800-7600 


TT Aqr 


2444131 


1979 Sep 15 


1.13 ± 0.02 


155 


4800-7600 


TT Aqr 


2444133 


1979 Sep 17 


1.16 ± 0.03 


156 


4800-7600 


TT Aqr 


2444138 


1979 Sept 22 


1.12 ± 0.02 


157 


4800-7600 


TT Aqr 


2444439 


1980 Jul 19 


1.31 ± 0.12 


159 


4800-8200 



Note. — Summary of polarization measurements of 60 Cyg and tt Aqr obtained between 1979-1987 using the Lyot Polarimeter at the 
University of Wisconsin's Pine Bluff Observatory 36" telescope. 
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Table 4 
Field Star Polarization Data 



Target Star 


Field Star 


Spectral Type 


Distance(pc) 


%Pol. 


P.A. 


TT Aqr 


HD213789 


G6III 


136 


0.53^ 


133.7^ 


TT Aqr 


HD213119 


K5III 


178 


0.23^ 


IIO.4I 


TT Aqr 


HD212320 


G6III 


141 


0.44^ 


I6O.5I 


TT Aqr 


HD211924 


B5IV 


282 


1.16^ 


140.21 


TT Aqr 


HD211838 


B8III 


228 


0.46^ 


138.8^ 


TT Aqr 


HD211304 


B9 


306 


0.221 


107.6^ 


TT Aqr 


HD211099 


B9 


455 


0.26^ 


119.2^ 


60 Cyg 


HD198915 


B6V 


288 


0.192 


15.02 



Note. — Note th at spectral types were obtained from SIMBAD, and distance measurements were obtained from the Hipparcos catalog 
llPerrvman et anil9 97). For reference, the Hipparcos der i ved d istanc es to n Aqr a nd 60 Cyg are 340pc and 418pc respectively. Archival 
polarization data were compiled from ^ IMathewson fc FbrdI 1IT970I ') and - jHeilej ll2000h . 
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Table 5 

Summary of Interstellar Polarization Values 



Star 


Value 




^max 


PA 


K 


C 


60 Cyg 


q' 


-0.076 ± 0.000 


6977 ± 241 


135.0 


1.17 





60 Cyg 


u' 


-0.083 ± 0.001 


6977 ± 241 


90.0 


1.17 





60 Cyg 


p 


0.112 ± 0.001 


6977 ± 241 


41.5 


1.17 





Pi Aqr 


q' 


-0.265 ± 0.001 


4959 ± 15 


135.0 


0.83 


0.04 


Pi Aqr 


u' 


-0.462 ± 0.001 


4959 ± 15 


90.0 


0.83 





Pi Aqr 


p 


0.514 ± 0.001 


4959 ± 15 


108.7 


0.83 






Note. — The derived modified Serkowski law parameters (jSerkowski et aLlll975l : IWilking et al.|[l982l ) derived for the perpendicular (U ) and 
parallel (Q ) ISP components. The total ISP estimate for each star was computed by vectoral addition of the quoted parallel and perpendicular 
components. Note that column 7, "C", denotes a constant term added to the modified Serkowski-law fit of the tt Aqr parallel ISP component 
to accomodate the presence of a possible small intrinsic polarization component. 
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Fig. 1. — The time evolution of 60 Cyg's observed (intrinsic and ISP) V-band polarization (top panels) and Hq EW (bottom panels) 
is illustrated for our entire dataset (left figure set) and a zoomed region (right figure set). Open circles correspond to Ritter observations 
while open triangles correspond to HPOL observations, and error bars (compiled in Table [T]l are smaller than the size of the data points. 
As discussed in Section |2] a factor of -1.9 A has been added to HPOL EW data compiled in Table [T] The three vertical lines in the right 
figure set depict (from left to right) the epoch of the maximum V-band polarization (JD=2451049), the maximum Hq EW (JD=2451171), 
and the minimum Ha EW (JD=2452038). The monotonically decreasing Ha EW strength demonstrates the gradual loss of 60 Cyg's disk 
over a time-scale of ~870 days. The ~120 day lag of the maximum Hq EW from the maximum V-band polarization and the ~190 day lag 
of the onset of the minimum Ha EW from the onset of the minimum V-band polarization suggest that the disk-clearing progressed in an 
"inside-out" manner. The full time-scale for disk loss, from the maximum V-band polarization to the minimum Ha EW, is ~1000 days (^870 
days -|- ~120 days). 
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Fig. 2. — The temporal evolution of tt Aquarii's observed (intrinsic and ISP) V-band polarization (top panels) and Ha EW (bottom panels) 
are illustrated for our entire dataset (left figure set) and a zoomed region (rig ht figure set). Open tr iangles correspond to our HPOL data while 
open circles correspond to literature V-band polarization values compiled bv lMcDavidI II1986I . [199 0. 1994, 1999). HPOL EW and polarization 
errors (compiled in Table [2]l are smaller than the size of the data points. The vertical lines in the right figure set denote (from left to right) 
the onset (JD=2448226) and termination (JD=2448612) of the ~390 day long first polarimetric outburst and the onset (JD=2448979) and 
termination (JD=2449713) of the ~730 day long second polarimetric outburst. 
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Fig. 3. — The polarization values of stars nearby the locat ion of 60 Cyg, based on the polarization catalog of lHeilej II2000IV One star having 
a similar distance (250 - 550 pc) as 60 Cyg (418 pc Perry man et al. ,: 199 7i[). shown with a solid line, exhibits significantly different polarization 
than 60 Cyg, indicating that a useful ISP estimate along the line of sight to 60 Cyg cannot be obtained from these data. For reference, the 
length of the vectors in the legend correspond to a polarization magnitude of 0.325%. 
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Fig. 4. — The V-band polarization of our entire dataset for 60 Cyg is plotted on a Stokes Q-U diagram. Note that the polarization errors 
(compiled in Table[TJ are smaller than the size of the data points. The data follow a strict linear trend indicating a constant scattering angle 
(i.e. polarization position angle) as expected in an electron scattering gaseous circumstellar disk system. The slope of the best fit linear 
regression to the data, 0.71 ± 0.02, yields the intrinsic polarization position angle of 60 Cyg, = 107.7 ± 0.4°. This implies that 60 Cyg's 
disk is oriented at a position angle of 9disft = 17.7° on the sky. 
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Fig. 5. — The V-band polarization of our entire dataset for tt Aqr is plotted on a Stokes Q-U diagram. Note that the polarization errors 
(compiled in Table [2]| are smaller than the size of the data points. Similar to 60 Cyg (Figure |4]l, the data follow a linear trend indicative of 
a constant scattering angle as expected in a circumstellar disk system. The slope of the best fit linear regression to the data, -0.50 it 0.01, 
yields the intrinsic polarization position angle of tt Aqr, 9» = 166.7 ± 0.1°. This implies that tt Aqr's disk is oriented at a position angle of 
Sdisk = 76.7° on the sky. 
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Fig. 6. — The wavelength dependence of 60 Cyg's polarization low-state data, binned to a constant (Poisson photon statistic) polarization 
error of 0.01% (depicted as a representative error bar), in the rotated Stokes U component (see Section 13. 3.31 1 was fit with a modified 
Serkowski law l lSerkowski et aLlll975l : rWilkinE et al.|[T982l) . yielding the estimate of the perpendicular ISP component along the line of sight to 

the star. The data in the rotated Stokes Q component, representing the parallel component of the ISP, was close to zero and not well fit by 
a Serkowski law. Section |3. 3. 3l provides a discussion of how we estimated the parallel ISP component for 60 Cyg. The parallel, perpendicular, 
and total ISP values adopted in this paper are summarized in Table [5] 
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Fig. 7. — The low state (diskless) 60 Cyg V-band polarimetric data are plotted on a Stokes Q -U diagram. Errors are smaller than the size 
of the data points. As discussed in Section l3.3.3l and illustrated in Figure|6] we fit a modified Serkowski law to the U component of these data 
to extract the perpendicular component of the interstellar polarization (denoted here by the solid line) . Since the wavelength dependence of 
the U data were not well represented by a Serkowski law (see Figure |6] we determined the minimum parallel component of the interstellar 
polarization to be simply the vector along the Q direction which joined the end of the perpendicular ISP component with the data having 
the minumum Q value, (denoted by the dashed line). This minimum parallel ISP component is tabulated in Table [S] 
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Fig. 8. — The residuals of our low-state 60 Cyg data, binned to a constant (Poisson photon statistic) polarization error of 0.015%, are 
shown after subtracting the total ISP values listed in Table [5] The data are consistent with Q = 0% and U = 0% and exhibit no significant 
wavelength dependence, providing confidence in the total interstellar polarization values we derived for 60 Cyg. 
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Fig. 9. — The wavelength dependence of tt Aqr's polarization low-state data, binned to a constant (Poisson photon statistic) polarization error 
of 0.01%, in the rotated Stokes U component (see Section |3. 3.31 1 was fit with a modified Serkowski law l lSerkowski et aLlll975l : FWilking et al.l 

119821 ). yielding the estimate of the perpendicular ISP component along the line of sight to the star (shown). The data in the rotated Stokes Q 
component, representing the parallel component of the ISP, were also fit with a modified Serkowski law (shown), which included an additional 
constant term to account for the presence of a potential remnant intrinsic polarization component. The parallel, perpendicular, and total ISP 
values adopted in this paper are summarized in Table [5] 
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Fig. 10. — The residuals of our polarization low-state n Aqr data, binned to a constant (Poisson photon statistic) polarization error of 
0.01%, are shown after subtracting the total ISP values listed in Table [5] The data exhibit a low level of residual polarization (< 0.05%) in 
both the Stokes Q and U parameters and exhibits little structure outside of the regions in which the HPOL CCD sensitivity is low (< ~4000 
A and > '^9500 A, providing confidence in the total interstellar polarization values we derived for n Aqr. 
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Fig. 11. — The absolute value of the error weighted deviation (deviance) of every V-band polarization observation of 60 Cyg (top panel) and 
TT Aqr (bottom panel) from the best fit line in each star's Q-U diagram (Figures|5]and|4{ is plotted as a function of the magnitude of intrinsic 
polarization present in each observation. Filled triangles correspond to observations during the low (diskless) state of each star, while open 
diamonds correspond to observations during the high (disk present) state. For both 60 Cyg and n Aqr, the likelihood of observing a small-scale 
deviation in the polarization PA (from the bulk intrinsic PA of the disk) clearly increases with the magnitude of intrinsic polarization present. 
A simple false discovery rate analysis (see Section |4.2|I confirms that this trend is statistically significant. 
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Fig. 12. — The absolute value of the error weighted deviation (deviance) of every V-band polarization observation of tt Aqr from the best 
fit line in the star's Q-U diagram (Figure |5ll is plotted as a function of time (bottom panels) along with the V-band polarization as a function 
of time (top panel). The right-side set of figures, a zoomed view of the entire (left-side figures) dataset, illustrates that the majority of the 
large deviations in the deviance correspond to one of the major polarimetric outbursts observed in n Aqr. The rise of the first component of 
the outburst (points A to B, < 54 days) is similar to that of the second component of the outburst (points D to E, < 35 days), while the 
time-scale of the decline phase in each component is also similar (points B to C, < 13 days; points E to F, < 10 days). We interpret these 
data as evidence of departures from axisymmetry in the disk and postulate that they may be evident of newly injected disk material running 
into a density structure such as a spiral densit y wave (Figure [141 . o r of the injection and subsequent circularization of new blobs into the 
inner disk, similar to the phenomenon noted bv lCarciofi et akl II2007I ). 
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Fig. 13. — The absolute value of the error weighted deviation (deviance) of every V-band polarization observation of 60 Cyg from the best 
fit line in the star's Q-U diagram (Figure [J]! is plotted as a function of time (bottom panel) along with the V-band polarization as a function 
of time (top panel). The right-side set of figures, a zoomed view of the entire (left-side figures) dataset, illustrates that the majority of the 
large deviations in the deviance correspond to large polarimetric outburst observed in 60 Cyg. Both the time-scale of the rise in the outburst 
(points A to B, < 26 days) and the time-scale of the decline in the outburst (points B to C, < 14 days) are similar to the time-scales in the 
analogous phases observed for tt Aqr. We interpret these data as evidence of departures from axisymmetry in the disk and postulate that 
they may evidence of newly injected disk material running into a density structure such as a spiral de nsity wave (Fig ure [T4ll or of the injection 
and subsequent circularization of new blobs into the inner disk, similar to the phenomenon noted bv lCarciofi et all teOOTD . 
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Fig. 14. — a) A simple cartoon image depicting au edgc-on Be star disk, whereby the inner region of the disk responsible for producing 
the observed intrinsic polarization is shaded yellow and the bulk of the region responsible for producing the observed Ho emission is shaded 
red. In the standard case of an axisymmetric disk, the maximum polarization vectors are oriented 90° from the plane of the geometrically 
thin disk, b) Departures from axisymmetry in the inner-region of the disk, such as those caused by the inject of one or more blobs which 
circularize into rings, would lead to intrinsic polarization vectors which are not all aligned perpendicularily with the bulk PA of the disk. We 
therefore speculate that one possible cause of the temporary (1-2 month) deviations in the intrinsic polarization PA values of both vr Aqr and 
60 Cyg is that blobs newly injected into the disk are created a warp-like structure in the inner disk, c) Alternatively, the observed deviations 
from axisymmetry could be caused by a blob (or series of blobs) launched from a non-equatorial latitude from the stellar photosphere, leading 
to the formation of a ring in an inclined (non-coplanar) orbit with respect to the semi-major axis of pre-existing disk material. 



